to/) 
3 



Heat Kernel Analysis for 
Ornstein-Uhlenbeck Operators with 
Quadratic Potentials 
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Abstract 

In this paper, we study Ornstein-Uhlenbeck operators with quadratic 
potentials. We use Hamiltonian formalism to characterise the singu- 
larities produced by the potentials by finding explicit geodesies of the 
operators, and obtain the heat kernels via a probabilistic ansatz. All 
the formulae are closed. 
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1 Introduction 

1.1 Background 

Ornstein-Uhlenbeck operator is a Linear operator 

T = -A + x- V 

acting on L 2 (M. n ; fi n ), where 

^ {E) = (2^72 // X P (-4) dX 

is the Gaussian measure, A = Y^j=i ®x an d V = (d Xl , d X2) . . . , d Xn ) denote 
Laplace operator and gradient, and x ■ y = YTj=\ x jVj * s the canonical inner 
product in R™ with an extension x ■ V = J2]=i x jd Xj - 
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For s > and / G L 2 (W l ; //„), 

e~ sT f(x)= [ f(e- s x + Vl-e^y)d fJ , n (y) 

is the Mehler formulae for the associated Ornstein-Uhlenbeck semigroup 
{ e_sT } s>0 - Readers may consult |Bog98], |FP94j and |Har05] for more el- 



ementary properties on Ornstein-Uhlenbeck operator and semigroup, and 
|Har99] . [F73] . [N73] for further interests. 

Ornstein-Uhlenbeck operators in the whole space L p (M. n ) are well under- 
stood (cf. [MOT], |MPV05] and |MPRSn2] ). and for operators on domains, 
bounded or unbounded, we refer |DPL04] . [FMP04| and [GHHW05] for fur- 
ther reading. Ornstein-Uhlenbeck semigroups also have rich applications to 
inequalities in geometry and analysis, such as Brascamp-Lieb inequalities (cf. 
[Bal89], |Bal91] and |BarC-E0"4] ). Gaussian logarithmic Sobolev inequality 
(cf. [L03]) and the reverse Bobkov isoperimetric inequality (cf. |BarC-EF01| 
and |Bob97] ). 



1.2 Motivation, methodology and arrangements 

In this paper, we shall consider Ornstein-Uhlenbeck operators perturbed by 
quadratic potentials. Essentially, they take the form 

L = —9d? + axd x + px 2 . (1.1) 

with all the parameters real scalars. The study of Hermite type operators (cf. 
|CF11] and |F12j ) show us that negative potentials give rise to singularities. 
We are naturally concerned with the following questions. 

• How do singularities arise? 

• How to characterize the singularities? 

Similar to Ornstein-Uhlenbeck operator, Kolmogorov operator in [K34J 
was first proposed to describe the probability density of a free system, and its 
heat kernel was obtained via Fourier transform. Hormander generalised the 
method to an appreciable extent in the introduction of his celebrated paper 
|Hor67j . where heat kernels of Ornstein-Uhlenbeck operators were benefited. 
However, operator in (11. ip is invariant under Fourier transform, which make 
those methods inapplicable to find its heat kernel. Here, we adopt Hamilto- 
nian formalism [BeaGG97j and an ansatz [Bca99j to reach our goals. 
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In section [21 we use Hamiltonian formalism to find the singularities and 
explicit geodesies of the operators ( II. ip . Heat kernels are obtained, with the 
help of ansatz, in section [3j We stress the critical case for singularities to 
arise and make more discussions on heat kernel in section HI 

As typical examples, perturbed Ornstein-Uhlenbeck operators 

L+ = -d 2 x + xd x + x 2 (1.2) 

and 

IT = -d 2 x + xd x - x 2 (1.3) 
will help readers digest our presentation. 



2 Hamiltonian formalism 



2.1 Hamiltonian systems 

For more generalities, we consider 1-d operators 

L = -6d 2 x + (ax + b)d x + px 2 , (2.1) 

where 6 > 0, a, b G R, p G K \ {0}. The Hamiltonian of the operator is 
defined as its full symbol 

H = -6? + (ax + b)i + px 2 , 

and the associated Hamiltonian system is 

dH _ 

dH 



dx 



ax -26£ + b 
= —2px — a£ 



Equivalent ly, 

Denoting 

we have 
with 



X 




a -26 




X 




~b 






—2p —a 






+ 






X 



x 



X = AX + B 



A 



a 

-2p 



-29 
-a 



B 



b 




(2.2) 



(2.3) 
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2.2 Homogeneous equation X = AX 

In order to solve the Hamiltonian system, we first consider the fundamental 
matrix of the system, i.e. solutions of the homogeneous equation 



X = AX. 
By the characteristic equation of A 

A - a 26 

\Xl-A\ = 

2p A + a 



(2.4) 



A 2 - a 2 - Ap6, 



we have two eigenvalues, Ai = \J a? + 4p6 and A 2 = — \/ a? + 4p6>. According 
to the signs of the quantity a 2 + 4/96*, namely positive, zero and negative, we 
find the fundamental matrix case by case. 



Case 1 a 2 + Ap6 > 0. (Ai > > A 2 ) 

Fundamental matrix $(s) has the form [e AlS ?7i, e A2S r] 2 ], where rji is an eigen- 
vector corresponding to the eigenvalue A«. To find 7/j's, we consider the 
characteristic equations 

(AJ-A)y = 



"Ai-a 26 




2/i 




"0" 


2p Aj + a 




2/2 








Taking 



26* 



Ai-a' 



26> 



2/1 



A,-a 
y 2 GK\{0} 



2/2 



, ^2 



26* 



A 2 -a> 



20_ e X lS 



Xi—a 



\is 



, we find a fundamental matrix 

(2.5) 



20_ e A 2 s- 



A 2 — a 
e X 2 s 
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The inverse of the fundamental matrix 

(Ai - a)(A 2 - a) 



20(X 1 - A 2 ) 



gA2S _26_gA2S 
A2— a 



_ g Ais __2#_gAis 
Ai— a 



(2.6) 



will be used later. 



P_ 

Ai 



g A 2 s 20 gA^s 

A2— a 



Ai— a 



(2.7) 



Case 2 a 2 + 4p0 = 0. (A x = A 2 = 0) 
Noting A is nilpotent of step 2, we may take 



Similarly, 



$(s) = exp(sA) = I + sA 



$ *(s) = exp(-sA) = 



1 + as -20s ' 
— 2ps 1 — as 



1 - as 20s ' 
2ps 1 + as 



(2.8) 



(2.9) 



Case 3 a 2 + 4p0 < 0. (Ai, A 2 e tR, 1 = ^1) 



In this case, Ai and A 2 are pure imaginaries. Denoting A = \f— a 2 — 4p6, 
we write Ai = i\o and Ai = —tXo- As same as case 1, we have fundamental 
matrix 



with 



and 



29 c \,s 26 c A 2 s' 

Ai— a, X'2~a 



D Xis 



*(0) 



29 



3 A2S 



2(9 



Ai— a A2— a 
1 1 



^(o) 



(Ai - a)(A 2 - a) 
20(Ai - A 2 ) 



1 



29 n 



A 2 — a 



26> 
Ai— a. 
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Given a real system (12 .4p . we expect to find a real fundamental matrix via 
the complex matrix ^f(s). We take 



$(s) = exp(sA) 

_ (Ai - a)(A 2 - a) 
26(X 1 - A 2 ) 



26 c Ais 26 e A 2 s' 

Ai— a A2— a 



3 Ai* 



= A2i 





" 1 


26 -1 

A2— a 




-1 


26 




Ai— a. 



cos(Aqs) + ^ sin(A s) -g sin(A s) 



-|£sin(A oS ) 



cos(A s) - ^ sin(A s 



(2.10) 



A direct computation gives 
$ -1 (s) = exp(-sA) 



26 sin(A .s) 



'cos(Aos) - % sin(A s) Aq 

|^sin(A s) cos(Aqs) + ^sin(Aos) 



(2.11) 



2.3 Geodesies 

Geodesies of the Hamiltonian system (12 ,2p is the solutions of 

x = ax + b q> 

with the boundary condition 

z(0) = x , x(1) = x. (2.12) 
The solutions of the equation (I2.3P has the form 

X{s) = $(s)C(s) (2.13) 
where $(s) is a fundamental matrix of the homogeneous equation 

x = ax. (J22D 
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and C(s) is a column vector satisfying the matrix equation 

Hs)C(s) = B, 

i.e. 

C(s) = ^- 1 (s)B. (2.14) 

We denote C(s) = [Ci(s), C 2 (s)]*, and c = [ci,c 2 ]' is a constant vector. We 
next identify the geodesies case by case. 



Case 1 a 2 + 4p9 > 0. (Ai > > A 2 ) 
By (F22D and fT2TT4|) . we have 
[Ci" 



C 2 



$ _1 (s)£ 



g A2S _2S_gA2i 
A2— a 



_ g Ais 26* g Ais 

Ai — a 



Hence, 



pb 
Ai 



D A2S 



Ci(«) 

c 2 (*) 



^e A2S + ci 



A1A2 



Making use of (I2.13P and (12. 5p . one has 



*(S)C{8) 



29 c Ms _J2e_ e \ 2 s 
Ai— a A2— a 



3 A2S 



^e A2S + a 



Ai A2 



ff e AlS + c 2 



._|i_ e A, ^A. + Ci _ _^ e A 2S | e A lS + 



e AlS ( -T4re A2S + ci ) + e A2S ( f|e AlS + c 2 



AiA 2 



7 



Thus, 



x(s) 



2bp9 



2bp9 



A X A 2 (Ai - a) A? (A 2 - a) 



3 (Ai+A 2 )s 



2 Cl 6 
X 1 — a 



2 C ,e 



A 2 - a 



2abp6 (Ax - A 2 ) A AaS 

~2— rr ryr r + «ie + a 2 e 

AfA 2 (Ai - a) (A 2 - a) 



(2.15) 



where rfi = — y 2 ^- and <i 2 



2c 2 6> 



Ai — a ^ A2— a' 

Putting A = \/a 2 + 4^^, we have 
Ai = A , A 2 = — A , 



2abpQ (Ai - A 2 ) 



ab 



AfA 2 (Ai - a) (A 2 - a) A§" 
The boundary conditions (I2.12p force rfj's to satisfy the linear equations 



" 1 1 ' 






e Ai e A 2 




d 2 



x + #' 



A o J 



One substitutes the solution 



di = 

d 2 p-Ao — pA 



-e^(xo + f|j + (x + f| 



into (I2.15P and finally obtains the explicit geodesies 



ab ( ab\ sinh (An(l — s)) ( ab\ sinh (Ans) 



AH ' V A? 



sinh Ac 



(2.16) 



Remark 2.1. It is apparent that p > ensures a 2 +4p# > 0, and hence there 
exists unique geodesies connecting two given points. This fact coincides well 
with Hermitian case, where positive potentials keep the geodesies regular. 
However, p < do not necessarily produce singularities in our case, as long 
as a 2 + 4p6 > 0. 



S 



Case 2 a 2 + 4p9 = 0. (A x = A 2 = 0) 
By fl23D and (12TT4]) . we have 

C = $- 1 (s)S 



~ci 




1 - as 29s 




~b 


A 




2ps 1 + as 








CM 

C2(S) 



bs - f s 2 + ci 
pbs 2 + c 2 



Hence, by (J2H3D and ([23]), one has 



1 + as -29s 
—2ps 1 — as 



'bs - f s 2 + ci 
pos 2 + c 2 



(1 + as) (&s - fs 2 + ci) - 2#s (p£>s 2 + c 2 ) ' 
-2ps (bs - fs 2 + ci) + (1 - as) (pbs 2 + c 2 ) 



Noting that a 2 + 4p9 = 0, one gets 



x(s) = —s + d 2 s + di 



(2.17) 



where d\ = c\ and d 2 = b + ac\ — 29c 2 satisfy the linear equations implied by 
the boundary conditions (12.121) 



d\ — xq 

ah 

d\ + a 2 = x - — 



d\ — Xq 

ab ■ 

d 2 = X - Xq - — 



After substituting dj's into f)2.17p . we finally obtain the geodesies 



ah 



x(s 



ab 



— S + \ X — Xq — 1 S + X , < S < 1. 



(2.18) 
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Case 3 a 2 + Ap9 < 0. (A 1; A 2 G iR, % 
By fl2TTD and we have 



1) 



c 2 



cos(Aqs) — ^- sin(Aos) 



f sin(Aos) 



2j sin(A s) cos(Aqs) + ^ sin(A s) 



Ci(s) 



^ sin (A s) + f| cos (A s) + ci 



■^cos(A s) + c 2 



Hence, by flgJSfl and (IZTOjl . one has 

= $(,)C( S ) 



cos(A s) + ^ sin(A s) sin(A s) 



-^sin(A oS ) 



cos(Aqs) - ^ sin(A s) 



6 sin (A s) + f| cos (A s) + c x ' 



^cos(A s) + c 2 



and an arrangement gives 



x(s) = ( cos(Aos) + — sin(Aos^ 
Ao 



b_ 

Ao 



sin (A s) + ^ cos (^o«) + ci 
A o 



— sm(A s) — cos (A Q s) + c 2 

Ao V A 



ab ( a . /x A 2#sin(A s) 
+ I cos(A s) + — sm(A s) I ci c 2 



A, 



An 



with Cj's fulfilling the boundary conditions f)2.12p 



1 







C\ 




x - 


ab 
A 


cos Ao + sin Aq 


28 sin Ao 
Ao J 




. C 2. 




x — 


ab 
\ 2 



(2.19) 



(2.20) 
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• sinA 7^ 0, {ci} i=1 have a unique solution. 
Indeed, 



1 







cosAo + ^sinA 



-1 




ab 




x - 


A 




x — 


ab 







ab 

XQ — T5 



Ao 



29 sin Aq 



(cos A + ^ sin A ) (x 



ab 



Ao 



29 sin Ao 



_ ab 
X -TJ 



Substituting q's back to (12 . 1 9[) and making an arrangement, one formulates 
the geodesies 



ab 



ab\ sin (A (l — s)) 



sin Ar 



+ \ x - 



ab\ sin (A s) 



An / sin A 



< s < 1. 
(2.21) 



• sinAo = 0, {q} 2 =1 have no solution or uncountably many solutions. 
Noting Aq = kir and the coefficient matrices of equations f)2.20p 



cos A + sin A 



one concludes that 







29 sin Aq 
A 



ab 

X 0~ T2 



_ ab 

X — JJ 
A 







-ir o 



aft 

^0 - 



A J 



=1 are solvable 
ab 



r ~ ^| = (~l) fc 

a b r -,\k 

> - ~ = (-1) I X 



A 



ab 



k 2 n 2 \ fc 2 7T 2 

x = xo, for & even 
2a6 



x 



k 2 n 2 



— xo, for k odd 
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and the solutions are 

ab 

Cl = x — 
_c 2 Gl 

By (12.191) . the geodesies are 



x ( s ) = + ( x o ~ J f cos(A s) + y sin(A s) ) - c 2 sin (A s) , (2.22) 



where c 2 is an arbitrary real scalar, and the parameter < s < 1. 

Remark 2.2. Condition sinAo = means that s = 1 is a singular regime, 
which corresponds to a hyperplane in the 2-d space-time coordinate (x, s). 
A detailed discussion for more singular cases has been made in |CFllj and 
|F12j . We do not pursue this point here. 

Finally, we summarise the explicit characterization of the geodesies as the 
following 

Theorem 2.1. Suppose that 9 > 0, a, b G 1R ; p G R \ {0}. The geodesies 
of the perturbed Ornstein-Uhlenbeck operators L = —Od^. + (ax + b)d x + px 2 
with the form {x(s)\0 < s < 1, x(0) = x , x(l) = x} are given by 



1. a 2 + 4p6 > 0. (a := V 'a 2 + 4pfl) 



ab f ab\ sinh (Aq(1 — s)) f ab\ sinh (Aqs) 



x(a) = -Ag + r° + A?J sinhAo + { X+ X 2 J sinh A 



V A / sum ao V A / sum ao 

(EHB 



2. a 2 + 4: P 6 = 0. 



ab 9 / ab\ . . 

z(s) = ys 2 + ( x-xo-y) s + x . fl2TT8|) 



3. a 2 + Ap6 < 0. (a := \/ -a 2 - Ap6 
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A ^kn (ke Z + ) 



ab f ab\ sin(A (l — s)) / ab\ sin (A s 



Aq V KJ shiXq \ Ag/ smA 

(E2U) 

A = kn (ke Z+) andx = ^ + (-l) fc (x - 

= + ( x o - £^ J ^cos(A s) + sin(A s)^ - c 2 sin (A s) 

flZSZD 

wii/i C2 an arbitrary real scalar. 



In particular, we obtain the geodesies of L 1 *" in (11.21) and (jl.3p . 

Corollary 2.1. There exist unique geodesies for the perturbed Ornstein- 
Uhlenbeck operator L + = —d 2 + xd x + x 2 fulfilling the boundary conditions 
x(0) = x , x(l) = x with 

x(g) = sinh(y^(l- g )) Xo+ sinh(V5 ! ) X) ^ 
sinh v 5 sinh v 5 



There exist unique geodesies for the perturbed Ornstein-Uhlenbeck operator 
L~ = —d 2 + xd x — x 2 fulfilling the boundary conditions x(0) = x , x(l) = x 
with 

. . sin(v3(l — s)) sin(\/3s) 
x(s) = — K y " x Q + — ^^x, < s < 1. (2.24) 
sin\/3 sin a/3 



3 Heat kernels 



3.1 A probabilistic ansatz 

We adopt a probabilistic ansatz to handle the heat kernel. It was first used 
in |Bea99] for Kolmogorov operators, and recently for operators of Hermite 
type |F12j . The perturbed Ornstein-Uhlenbeck operators 

L = -9d 2 x + (ax + b)d x + px 2 
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have heat kernels of the form 

P(t; xq, x) = ip(t) exp [a(t)x 2 + (3(t)xxo + 7(£)xq + p(t)x + z/(t)x } (3.1) 

where the coefficients tp, a, (3, 7, p, v are functions of the time variable t and 
satisfy the following systems of ordinary differential equations 



3.2 Identification of coefficients 

We now solve the equations fl3.2H - fl3.7j) and identify the free constant in each 
equation by comparing the coefficient with that in Hermitian case |F12j . i.e. 
a = 0. Without confusion, we denote the free constant by C which may be 
different from time to time. 

Case 1 a 2 + 4p9 > 0. 



Equation (13. 2ft : Putting a.\ = a+ ie + P and a 2 = - ^ P , one has 



a = 46a 2 — 2aa — p 

= 46afi - a/3 

1 = ep 2 

fi = 49ap — ap — 2ba 

v = 29 pp - b/3 

tp^tp = Op 2 + 29a - bp. 



(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 



a = 49(a — ai)(a — a 2 ) 




a — «2 




a 



v/a 2 + 4p9 1 + Ce 2 ^ a2+ipe 



Let (7=1, one gets 




(3.8) 
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Equation (13. 3ft : 

/T 1 /? = -yjc? + 4p#coth (V^ 2 + 4p6>t) 



ln|/?| = - In (sinh (^/a 2 +4p§tj] +C 
C 



P 



sinh [\fa 2 + 4p6»f 



Let C = ^ a ^ 4p9 , one gets 



20 sinh y/a 2 + 4p^t 



Equation (ETll : 



7 



a 2 + 



4^ sinh 2 ( y/a 2 + 4pfl t 



Let C = 4|, one gets 



a 2 + Ap6 



coth(Va 2 + 4p#t) +C. 



7(t) = £ - ^±^coth (v^T^t) . (3.10) 
Equation (13.51) : \i has the form 



sinh ^a 2 + 4p8 tj 
where C(t) satisfies 

C 



sinh (^/a 2 + Ap8 tj 



2ba(t). 
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By (I3.8p . an integration yields 

cosh ( ^a? + Ap9 sinh ( ^a? + Ap9 1 



C{t) = -2b 
Hence, 



AO^/a? + Ap9 4:0 



Equation (13. 6p : 

a ^ cosh (y/ a 2 + 4p6't 



2,9 sinh 2 (^a 2 + 4p0^ 
20 v /a 2 + 4p6>sinh i^a 2 + Ap9t 



+ C. 



Let C = one gets 



a6 6 
20 ^a 2 + 4p# sinh ( ^a 2 + Ap6 1) 16 



Equation (13. 7ft : By (13. 8 p and (13. lip , an arrangement of right hand side 
of (EH) gets 



V= -7- \ ~ 9 CQth (Va 2 + 4p0t 

Ad (a 2 + Apd) sinh 2 (y a 2 + 4p#t) 2 v 

(a Pb 2 \ 
\2 a 2 + A P e)' 



+ 
Hence, 



C 



sinh 1/2 L/a 2 +Ap9t 

exp { (f - ^T^) ' " coth (v 7 ^^) } 
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1/2 

) 



Let C = ( : ^Sr^l > one S ets 

1/2 



Va 2 + 4p6» 



47r#sinh L/a? + Ap9t 

6XP {(I ~ ^T^)'-^ + 4p6l)3/ 2 COth (v / « 2 + 4p^)}- 

(3.13) 



Case 2 a 2 + = 0. 
Equation (I3.2p : 



Let C = 0, one gets 



Equation (13. 3ft : 



Let C = t^, one gets 



Equation fl3T4j) : 



a — — 



-4#t + C. 



"4 



«*) = W (3-15) 



1 

7 



40 1 2 
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Let C = jjj, one gets 



. . a 1 
7(*) = T-o 



46> 4.9 1 
Equation f)3.5p : \i has the form 

/ \ C(t) 

where C{t) satisfies 

C 

— = -2ba(t). 
By (13.141) . an integration yields 



Hence, 



Equation (I3.6P : 



Let C — we, one gets 



s ab 9 6 



" W = "is* + 2*- 



ab 

v = ~Te 



^ , 



i/(t) = i H . 

W 49 29 



Equation ffXTj) : 



i a 2 o 2 9 1 /a b' 

(f ~V = — -r h 

^ ^ 169 2* V2 4i 



C f/a o 2 \ a 
V 3 = ~T. ex P ^ o _ 77i * 



1V2 40,/" ' 4 

1 /2 

Let C = life) ' one § ets 



1 \ 1/2 r / z ' 
1 \ ' f/a tr 
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Case 3 a 2 + 4p9 < 0. 
Put «o = \/~ a 2 — 4p6*. 
Equation 1ET2]) : 



■ / a \ 2 a o 



a a / „s 
a= A0 + ^ tan («o^ + C'). 



Let (7 = 1, one gets 



a V -a 2 - 4p6> 



= 40 ~ itf - ^ Cot ( V -« 2 - 4p0 * J • (3- 2 °) 



Equation H3.3|> : 

/? /3 = — a cot(ao^) 



c 



sin(ao t) 



Let C = H , one gets 



/?(*) = V- " 2 - 4 ^ r . (3.21; 



20 sin (v 7 -^ 2 -4p0t 

Equation (j37D : 



7 = 40 esc (a 1) 



7 = _ ^ cot («o^) + C- 



Let C = ^, one gets 



a v — a 2 — 4p# 



7 ^ = 4^ ~ Cot ( V -« 2 - 4p0 * J • (3-22) 
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Equation (I3.5p : \x has the form 

C(t) 



n(t) 



sin(ao t) ' 
where C(t) satisfies 

C 

-2ba(t). 



sin(a 1) 
By ( I3.20p . an integration yields 



„. . ab , . b 

Hence, 

u(t) = ab cot (j-a 2 -Ap9t) + — . 

Equation (13. 6p : 

a& cos(cko t) 



v 



26>sin 2 (a *) 



v - a ^ + C 

29a sin(a 1) 



Let C — A, one gets 



20 y/- a 2 - Ap9 sin f ^-a 2 - 4p6>i) 2y 



Equation (13. 7p : 



i a 2 6 2 2 a { a pb 2 

^ = « CSC (ao*)- T cot(a t)+ - + 



C (fa pb 2 \ a 2 b 2 , / 

* = exp 1 U + 7? J * " 4^I cot(aot) . 
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sin f^-a 2 -4p0*J y 

• 6XP { (f " ^fv) * " 49(-J-4p9)W COt (V-a 2 -4p(?t)}. 

(3.25) 



3.3 Closed heat kernels 



Making use of ansatz (13. ip . we summarise the results on coefficients in the 
following theorem. 

Theorem 3.1. Suppose that 9 > 0, a, b E R, p E R \ {0}. The heat kernels 
P(t; x,Xo) of the perturbed Ornstein-Uhlenbeck operators L = —9d 2 + [ax + 
b)d x + px 2 are given by 

1. a 2 + 4p8 > 0. 



P(t;x,x Q ) 



^a 2 + Ap9 



1/2 



47r6»sinh (^Ja 2 + Ap6t 
pb 2 



exp 



• exp 



exp 



exp 



exp 



2^2 



t - 



a 2 b 



2 a 2 + 4p6 J 46{a 2 + ApOf' 2 



coth (ya? + 4p6t 



a 2 + 4p6 



29 sinh (y« 2 + Ap9t 
ab 



XXq 



IB J a 2 + Ap9 



coth (^H4^i) + 



x 



ab 



26 y/ a 2 + 4p#sinh (^a 2 + 4p0t 



(3.26) 
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2. a 2 + Ap9 = 0. 




3. a 2 + Ap9<0. 



P(t;x,x ) 

^-a 2 - Ap9 
An9 sm (y '-a 2 - Ip6t) 



1/2 



cxp 



• cxp 



exp 



cxp 



• exp 



a 2 b 2 



\ - Jhie) 1 - AQ{-a 2 -Apeyi 2 cot < v ~ a2 ~ Apet 



^-a 2 - Ap9 
2e sin (^f-a 2 -ApOtj 

( , ^ cot (f^i) + *) 

\29^-a 2 -Ap9 V v J 29 J 



x 2 + x 2 ) 



x 



ab b 
29 ^ -a 2 -Ap9 sin {^/-a 2 -Ap9^j 26 



(3.28) 
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4 Further discussions on geodesies and heat 
kernels 



4.1 Critical cases 

The critical case, i.e. a 2 + Ap9 = is a watershed in the discussion of singu- 
larities. Due to its independent interest, we formulate it as a proposition. 

Proposition 4.1. Suppose that 9 > 0, a, b 6 R, p e R \ {0} and a 2 + ApO = 

0. Perturbed Ornstein-Uhlenbeck operators L = —9d 2 + (ax + b)d x + px 2 
have unique geodesies {x(s)\0 < s < 1} fulfilling boundary conditions x(0) = 
x , x(l) = x with 

oft , / ab\ . . 

x{s) = —s 2 + Ix -x - y) s + x o> 

and the heat kernels are given by 
P(t; x,x ) 



1 V /2 [fa b 2 \ a 2 b 2 



exp { \Te ~ Wt) ^ + xl) + Wt XXo + 



ab i. 

T9 t + Te ]{x + Xo < 



(E2ZD 



Remark 4.1. If ab ^ 0, there exists a unique parabola connecting x and x. 
If ab = 0, we may assume a ^ and 6 = since other cases are of Hermite 
type (cf. [F12] ). The operators in the family 

[L = -6d 2 x + axd x + px 2 \a 2 + 4p6 = 0} 

have exactly the same geodesic — a straight line connecting xq and x. In this 

2 

case, one infers that p = — jz < is the minimum value that p could take 
to keep the operators regular. In particular, perturbed Ornstein-Uhlenbeck 
operators L = —0d 2 + axd x + px 2 share the same geodesies with Laplace 
operator L = —d 2 , but the heat kernel of L is extinct from Gaussian. In 
other words, geodesies do not suffice to characterize heat kernels in our cases. 
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4.2 Heat kernels in higher dimensions 



For brevity, we assort the heat kernels in Theorem 13. 11 

Definition 4.1. For L = —9d 2 + (ax + b)d x + px 2 , formulae (13.26P is called 
heat kernels of Type I and denoted by Pi(t; x, xq). Similarly, (I3.27P and ( I3.28P 
are called Type II and Type III, and denoted by Pu(t; x, xq) and Pm(t; x, xq) 
respectively. 

We next formulate heat kernels for operators in higher dimensions, and 
illustrate them to L in higher dimensions. 

Proposition 4.2. Let 

n n n 

L Yl °i d xj + Z)( W + b i) d *i + Yl l'r'1 

3=1 3=1 3=1 

be a perturbed Ornstein-Uhlenbeck operator in n spatial variables, where 9j > 
0, dj, bj e K and pj el\0. The heat kernel P(t; x, Xq) is given by 

P(t; x, x ) = Piit] x, xo)P H (t; x, x )Pni(t; x, x ) 

Proof. Due to the independence of each spatial variable, the heat kernel P 
is a product of n 1-d heat kernels of the form f)3.26p . (I3.27P and (13.281) . 
According to the sign of the quantity a 2 + 4p6 and Definition 14.11 these heat 
kernels can be assorted into 3 parts, i.e. Pi, Pu and Pm, which completes 
the proof. □ 



Corollary 4.1. The n- dimensional perturbed Ornstein-Uhlenbeck operator 
L + = —A + x ■ V + Ixl 2 has heat kernel 



P(t;x,x ) 



n 

47rsinh(v / 5t) J 

{1 — \/5coth(\/5£) , 2 2N y/E 
4 (N +|l ° l) + 2smh(V5 t ) a:,I ° 

The n- dimensional perturbed Ornstein-Uhlenbeck operator L~ = —A + x 
V — \x\ 2 has heat kernel 

{1 — \/3cot(\/3t) /, ,9 .ON 
1 (M + M ' + ^Mvm x ' 10 
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5 Conclusions 



We first give direct answers to the questions proposed at the beginning of 
this paper. 

• The sign of the quantity a 2 + 4p8 decides whether singularities arise or 
not. Specifically, perturbed Ornstein-Uhlenbeck operators L = —9d 2 + 
axd x + px 2 is regular if a 2 + 4p9 > 0, and irregular if a 2 + Ap6 < 0. 
Remark 12.11 Remark 14.11 and Remark 12.21 make brief comparisons to 
Hermitian case on singularities. 

• We characterize the singularities from the point of geodesies and heat 
kernels. Complete discussions for each case are summarised in Theorem 
12.11 and Theorem 13.11 

An obvious extension of this result would be to obtain geodesies and heat 
kernels for the operators like 

L = -Qd 2 x + (ax + b)d x + p x 2 + p x x + p 2 . 

For operators in a more general form 

L = -(AV, V) + (Bx, V) + (Cx, x), 

this paper treats the diagonal case of the coefficient matrices. When A, B 
and C are not diagonal, some results have been reached. |Bea99] essentially 
obtain the results for A > and C > 0. However, C < is a specially 
interesting case to study, whose geometric analysis is absent in the literature. 
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